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Abstract — In this article we show the relation between the 
theory of pulse shaping for WSSUS channels and the notion 
of approximate eigenstructure for time-varying channels. We 
consider pulse shaping for a general signaling scheme, called 
Weyl-Heisenberg signaling, which includes OFDM with cyclic 
prefix and OFDM/OQAM. The pulse design problem in the 
view of optimal WSSUS-averaged SINR is an interplay between 
localization and "orthogonality". The localization problem itself 
can be expressed in terms of eigenvalues of localization operators 
and is intimately connected to the concept of approximate eigen- 
structure of LTV channel operators. In fact, on the £2-level both 
are equivalent as we will show. The concept of "orthogonality" in 
turn can be related to notion of tight frames. The right balance 
between these two sides is still an open problem. However, several 
statements on achievable values of certain localization measures 
and fundamental limits on SINR can already be made as will be 
shown in the paper. 

I. Introduction 

Optimal signaling through linear time-varying (LTV) chan- 
nels is a challenging task for future communication systems. 
For a particular realization of the time-varying channel opera- 
tor the transmitter and receiver design which avoids crosstalk 
between different time-frequency slots is related to "eigen- 
signaling" which simplifies much the information theoretic 
treatment of communication in dispersive channels. But it 
is well-known that for an ensemble of channels which are 
dispersive in time and frequency (doubly-dispersive) such 
a joint separation of the subchannels can not be achieved 
because the eigen decompositions can differ from one to 
another channel realization. Several approaches are proposed 
to describe "eigen'-signaling in some approximate sense. 
However, then a necessary prerequisite is the characterization 
of approximation errors. 

A typical scenario, present for example in wireless com- 
munication, is signaling through a random doubly-dispersive 
channel, represented by the operator T~L. The received signal 
r : R — > C at time instant t is then: 

r(i) = (Hs)(t) + n(t) 

A preferred design of the transmit signal s : R — * C needs 
knowledge on the true eigenstructure of the channel operator 
Ti. This would in principle allow interference-free transmis- 
sion and simple recovering algorithms of the information from 
received signal r degraded by the noise process n. However, 
for 7t being a random operator, a random eigenstructure has 



to be expected in general such that a joint design of the 
transmitter and the receiver for an ensemble of channels has 
to be performed. Nevertheless, interference can then not be 
avoided and remains in the communication chain. For such 
interference scenarios bounds on the distortion of a particular 
selected signaling scheme are mandatory for reliable system 
design. 

This problem has been considered, for example, already in 
[1] and [2]. The investigation of the £ p -norm E p of the error 
H.s — \r for A 6 C enables us to improve and generalize 
recent results in this direction. We will focus in this article 
on the results for p = 2 which show the important relation 
to pulse shape optimization for Weyl-Heisenberg signaling 
[3], [4], [5]. We will use this generalized description of a 
multicarrier system and consider in this framework the channel 
averaged "signal to interference and noise ratio". To perform 
this average we will use the wide-sense stationary uncorrelated 
scattering (WSSUS) channel model. 

The paper is organized as follows. In Section II we will 
review the basics from time-frequency analysis, introduce the 
spreading representation of doubly-dispersive channels and 
consider the problem of controlling the approximation error 
E2. In Section III we will focus on the description of multi- 
carrier transmission and WSSUS pulse shape optimization. We 
demonstrate the relation to i?2- In Section IV we give several 
relations to the "size" (Lebesgue measure) \U\ of the support 
U in the time-frequency plane of the channel's spreading 
function. In the last part we will verify our framework with 
some numerical tests. 

II. Time- Varying Channels 

A. Time-Frequency Shifts 

The fundamental operations on a signal 7 : R — > C caused 
by the time-varying nature of typical wireless channels are 
time-frequency shifts. Let us denote a displacement in time 
by /xi and in frequency by /12 by the operator defined as: 

(S M7 )W:=e^ 2t 7 (t-Mi) (D 

(i is the imaginary unit and /i = (/Xi,/^)). These operators 
act isometrically on the £ p (R) spaces 1 , hence unitary on the 
Hilbert space /^(R)- The abstract action of onto a function 

1 £ p (K n ) denotes Lebesgue spaces on M 71 with usual norms ||-|| p for 
1 < p < 00 



7 is observable for example as a correlation response with 
another function g. In this way the (cross) ambiguity function: 

A ff » d ^ f (g, S^) = f g(xh(x - ^e^'dx (2) 

JR 

can be defined as an observable quantity. Definition (1) is 
not unique and corresponds to a particular polarization in the 
Weyl-Heisenberg group (see [6], [7]). 

A fundamental operation in time-frequency analysis is the 
symplectic Fourier transform T a defined for a function / : 
R 2 -»■ C as 

(J>/)0*) = / e- i2 ^^f(u)du (3) 

JR 2 

where rj(u, p) := v\[i2 — V2[i\ is called the symplectic 
form. (3) differs from ordinary Fourier transform only by 
coordinate/sign switching which quite intuitive because time- 
shifts iii cause oscillations in frequency where frequency shifts 
\i2 cause them in time. Then the (cross) Wigner distribution 
is the symplectic Fourier transform of the (cross) ambiguity 
function. 

B. Spreading Representation of Doubly-Dispersive Channels 

Doubly-dispersive channels are physically characterized by 
scattering in the time-varying environment. Typically one 
defines therefore the channel operator Ti weakly in terms 
of the transmit and receive filters 7 and g of a particular 
communication device as: 

(ff,«7>=/ S„( M )A S7 ( M )d/i (4) 

JR 2 

The function £7^ : M 2 — > C is then called the spreading 
function of the channel operator H. Its symplectic Fourier 
transform, L-j-c = JF S S-^, is called the symbol of Ti. The 
practical assumption that the support U C M 2 of is 
compact renders Ti to be a Hilbert-Schmidt operator whenever 
S-h is bounded and \U\ > (single- or non-dispersive 
channels have \U\ = 0). In this paper we will restrict the 
analysis to those channels to avoid the use of distributions. 
To this end, let us call OP(Z7) as the set of channel operators 
TL having a spreading function Y^-j-c with non-zero and finite 
support in U where < |Z7| < 00. 

C. Approximate Eigenstructure 

It is a fundamental question, how well we can describe the 
action of the channel in terms of its symbol. More generally 
speaking, how valid is an approximation of the form: 

HS^j w A(/i)S , /i fi( (5) 

See for example [8] for recent applications of this approx- 
imation. For example, for g and 7 being singular functions 
of Ti (eigenfunctions of Tili.* and Ti.*7i) there is equality 
for a particular A(0) (the singular value). However, this is 
not well suited for our formulation, because Ti. will be a 
random channel operator as for example further described in 
Section II-D. However, if we restrict the set of operators to 
be from class OP(U) we can seek for g and 7 which are 



independent of the channel realization and reasonable bound 
the approximation error E p = \\TiS ^ — X(fi)S p g\\ p . For a 
single Ti we will call this setup as approximate eigenstructure 1 
and in this article we will consider p = 2 because of its relation 
to pulse design (considered in Section III). For more general 
results see [7]. Our problem is: 

Problem 1 Let be Ti E OP(U) and \\g\\ 2 = H7H2 = 1 and 
1 < a < 00. Let 5 E R + be independent of S^' such that 

E 2 < 5 ■ \\X H \\ a (6) 

How small can we choose 6 given g, 7, U and a ? What can 
be said about inf gn ($) given U and a. 

Note that this formulation comprises a joint approximate 
eigenstructure for all channels Ti G OP(f/) but with individual 
bounds 5 ■ ||E^|| a for the approximation error. 

First decompositions results in this field can be found 
already in the literature on pseudodifferential operators [9], 
[6]. More recent results with direct application to time-varying 
channels were obtained by Kozek [3] and Matz [10] which 
resemble the notion of underspread channels. They found 
results on 8 for A = L-j^ and a = 1 which follow from the 
approximate product rule in terms of Weyl symbols. This es- 
timate intimately scales with \U\ and breaks down at a certain 
critical size. Channels below this critical size are called in their 
terminology underspread, otherwise overspread. However, the 
critical term in their estimate can be dropped and the overall 
result can be improved to S 2 < 2||(1 — A gl )xu\\oo [7]. Similar 
estimates can be obtained for other weight functions instead 
of xu ( me characteristic function of U). For this case (a = 1) 
it is difficult to separate explicit relations on \U\. However, as 
will be shown, for a > 1 this can be done. 

With the choice A = L-h a formulation independent of g 
and 7 is achieved. An approach which is better suited for our 
application instead is 

A(/x) = (S^g, H, S^) = f s (E„A ST )M (7) 

This choice is also known as the orthogonal distortion 
(smoothing with the cross Wigner function) and corresponds 
to the E2 minimizer in (6) for fixed g and 7. This approach 
was also already considered for g — 7 in [2]. One can prove 
the following Lemma [7]: 

Lemma 1 Let 1 < a < 00 and 1/a + 1/6 = 1 and E-^ E 
£ Q (R 2 ). // \U\ < 1 it holds 

p / r- \ 1/ max(6,2) 

P^< \\U\- ^|A 9 »| 2 ^J (8) 

Similar bounds can be found for all p-norms with 2 < p < 00 
and more detailed results avoid also the constraint \U\ < 1. 

2 This approach is similar to the concept of a pseudospectrum. The set of 
numbers A for which exists a normalized 7 such that WHy — A7H2 < e is 
called the e-pseudospectrum of H. 



D. WSSUS Scattering 

A common statistical model for random doubly-dispersive 
channels is the WSSUS channel model, originally introduced 
by Bello [11]. In this model the spreading function S-^ is 
described as a zero-mean 2D-random process uncorrected at 
different arguments, which gives: 

E{|(. 9 ,-H 7 )| 2 }= / C( M )|A g »| 2 ^ (9) 

JR 2 

The (non-negative) function C : M 2 — > R + is called the 
scattering function. Let be ||C||i = 1. Then C(fi) is also 
the probability density for the occurrence of a time-frequency 
shift fj, in the channel. In time-frequency analysis the term 
(C, |A g7 | 2 ) in (9) is sometimes also called the modulation 2- 
norm [12] of 7 with respect to window q and weight function 

Vc. 

III. Signaling and Pulse Shaping 

In Lemma 1 and (9) we have already shown certain local- 
ization measures characterizing the signal distortion in time- 
varying channels, in particular for the WSSUS model. In this 
section we will now establish a simple communication chain, 
including OFDM and OQAM/OFDM, which use g and 7 as 
corresponding transmit and receiver filters. Finally, we aim at 
a relation between {g, 7}, the WSSUS-averaged SINR and the 
role of time-frequency localization which is relevant for £3. 

A. Weyl-Heisenberg Signaling 

Conventional OFDM and pulse shaped OQAM/OFDM can 
be formulated jointly within the concept of Weyl-Heisenberg 
signaling [13]. To avoid cumbersome notation we will adopt a 
two-dimensional index notation n = (ni, ri2) £ Z 2 for time- 
frequency slots n, such that the baseband transmit signal is 

S W = X! Xn1n{t) = X! X n (SAn 7 )(i) (10) 

where AZ 2 is a lattice (A denotes its 2 x 2 real generator 
matrix) with |A| := | dct(A)| > 0. The indices n range over 
the doubly-countable set I C Z 2 , referring to the data burst to 
be transmitted. In practice A is often restricted to be diagonal, 
i.e. A = diag(T, F). Other lattices are considered for example 
in [14]. The efficiency (in complex symbols) of the signaling 
is | TV | 1 . The coefficients x n are the complex data symbols 
at time instant m and subcarrier index ri2 (from now on 7 
always denotes complex conjugate and •* means conjugate 
transpose), where x = (. . . , x n , . . . ) T . By introducing the 
elements H mn :— (<? m ,7i 7 „) of the channel matrix H 
the multicarrier transmission can be formulated as the linear 
equation x = Hx + n, where h = (. . . , (g m , n), . . . ) T is 
the vector of the projected noise having variance a 2 := 
E{ | (.9m, n) | 2 } per component. 

B. Complex Signaling and OFDM 

The classical OFDM system exploiting a cyclic prefix 
(cp-OFDM) is obtained by assuming a lattice generated by 

A = diag(T,F) and ~ X[-T cp ,T u ](t) which is (up 



to normalization) the characteristic function of the interval 
[—T cp , T u ]. T u usually denotes the duration of the useful part 
of the signal and T cp is the length of the cyclic prefix such that 
T = T U + T cp . The OFDM subcarrier spacing is F = l/T u . 
At the receiver the rectangular pulse g(t) ~ X[o,T u ](t) is used 
which removes the cyclic prefix. The efficiency is given as 
| TV | 1 = T U /(T U + T cp ) < 1. It can be easily verified that 
A37 (/•* + Am) ~ 5 m> o if < fjhx < T cp and ^2=0 (or see 
[15] for the full formula). Thus, orthogonality is preserved 
in time-invariant channels with maximum delay smaller than 
T cp , but at the cost of signal power (the redundancy is not 
used) and efficiency. 

C. Real Signaling and OQAM 

For those schemes an inner product Re{(-, )} and A = 
diag(T, F) having |A| = 1/2 is considered, which is realized 
by OQAM based modulation for OFDM (OQAM/OFDM) 
[16]. Before modulation the mapping x n = i n x\ has to 
be applied 3 , where xf t G K is the real-valued information 
to transmit. After demodulation xf n = Re{i~ m x m } is per- 
formed. Moreover, the pulses (5,7) have to be real. Thus, 
formally the transmission of the real information vector a; R = 
(. . . , x^, . . . ) T can be written as x R = H R x R + n R where the 
real channel matrix elements are: 

Hf nn = Re{^™- m i^„ M^ } = Re{i n ~ m {g m ,H«( n )} (11) 

and "real-part" noise components are = Re{i~ m (g m , n)}. 
There exists no such orthogonality property for OQAM based 
signaling as for cp-OFDM in time-invariant channels, but 
biorthogonality of the form dl{(g m ,j n )} ~ & m ,n can be 
achieved. It is known that the design of (bi-)orthogonal 
OQAM signaling is implicitly related to (bi-) orthogonal Wil- 
son bases [17] which is in turn equivalent to the computation 
of (dual-) tight frames [18]. 

While the system operates with real information at e = 2 
the effective efficiency is one, which has advantages in the 
view of pulse shaping as will explained later on. 

D. A Lower Bound on the WSSUS Averaged SINR 

In multicarrier transmission most commonly one-tap equal- 
ization per time-frequency slot is considered, hence it is 
naturally to require E{|i7 mm | 2 } to be maximal and the 
averaged interference power from all other lattice points to 
be minimal as possible. For real schemes H m ^ m has to be 
replacedbyfl* n ftom(ll).For||C||i = 1, ||. 9 || 2 = H7II2 = 1 
and ~E{x m x n } = S m , n we have shown in [4] that the averaged 
"signal to interference and noise ratio" is lower bounded by: 

SINRfo t A]> ( C i l A ff7l 2 ) a2 ) 

with equality for {j m } forming a tight frame (more details and 
references on frame theory in [4]). The constant B 1 = p(5 7j a) 

3 We use here the notation i n = i n i+ n 2 Furthermore other phase 
mappings are possible, like j"i+™2+2n 1 n 2 



is the spectral radius of the positive semidefinite operator 5 7i a 
defined as: 

(S 7 , A /)(<):= J2 (SxlJ)(S xl )(t) 

AeAZ 2 

The minimal B 7 is max(|A| _1 , 1) which is achieved for tight 
frames in the case of |A| < 1 and for incomplete orthogonal 
bases in the case of |A| > 1. In addition to complex signaling 
as considered in Section III-B Equation (12) holds for real 
signaling (Section III-C) as well if the spreading function and 
the noise are circular-symmetric processes [4]. 

E. SINR Optimization and the Pulse Shaping Problem 

It follows that maximizing the bound in (12) is an interplay 
between two criteria: localization (maximizing (C, |A g7 | 2 )) 
and "orthogonality", respectively "tightness" (minimizing 
By). Orthogonalization is only possible for |A| > 1. For 
|A| < 1 a tight frame has to be constructed. However, with 
the notion of the adjoint lattice A° := |A| _1 A both cases are 
dual in the sense of Ron and Shen [19], which means that we 
can limit ourself to the case |A| < 1. 

One possible method [14], [4] to achieve a minimal B 7 
in a /^-optimal sense is to apply O := (|A| • S^a)" 2 on 
the result of a localization procedure 4 to obtain 7 = O7 
(for |A| > la similar procedure follows from replacing 
A by A°). For example, in [20] it was shown that IOTA 
(Isotropic Orthogonal Transform Algorithm, A = diag(T, F) 
with |A| = 1/2) as proposed in [21] is equivalent to O applied 
on a Gaussian. In combination with OQAM modulation as 
introduced in section III-C this is also an orthogonal signaling 
(in the absence of a channel). The advantage of this approach 
is that at |A| = 1/2 the operation O is quite stable. In contrast: 
optimal efficiency with orthogonal (or bi-orthogonal) complex 
signaling (this implies |A| = 1) is always affected by the 
Balian-Low Theorem (BLT) (see for example [22]), i.e. O is 
unbounded if it exists. 

Estimates on the localization penalty of the operation O 
on the SINR are in general related to mapping properties of 
O on modulation spaces [12]. Beside the general problem of 
invertibility of the operator in (III-D) (for recent results see 
here for example [23] or [22]) there remains the open question 
on optimal constants which are necessary to draw conclusions 
on the SINR. 

IV. Localization and Support Bounds 

A. Pulse Shaping and Approximate Eigenstructure 

The detailed knowledge of the scattering function of an 
statistical ensemble of channels can improve the overall per- 
formance of a communication device [4]. However, also this 
statistical parameters change overtime in realistic applications. 
A more reliable assumption is the support knowledge only, 
i.e. a scattering function which is the normalized indicator 
function C = Xu/\U\ of U. We then see that approximate 
eigenstructure in the formulation given in Lemma 1 and the 

_ 1 

4 S -r£ exists if {7n}„ e z2 is a frame. 



localization aspect (the term (C, |A g7 | 2 )) of pulse shaping 
(12) are the same problems. 

B. Localization Operators 

Since |A g7 | 2 is quadratic in 7 we can rewrite (|A ff7 | 2 , C) = 
(7, Lc.g^i) where this quadratic form defines (weakly) an 
operator Lc, g - Such operators are also called localization 
operators [24]. Lc, g is compact which follows from our 
assumptions on C, such that 

max(|A 97 | 2 , C) = A max (i c ,g) (13) 

7 

This means that for fixed g the optimal transmit pulse 7 
is given by an eigenfunction corresponding to the maximal 
eigenvalue of Lc, g - The latter can also be reversed, i.e. for 
fixed 7 the optimal receive filter g is determined. 

max(|A S7 | 2 , C) = A max (L 6>7 ) (14) 

where C(fi) = C(—fi). The eigenvalues and eigenfunctions 
of Gaussian (g is set to be a Gaussian) localization operators 
on the disc (U is a disc) or more generally with C hav- 
ing elliptical symmetry are known to be Hermite functions 
[24]. Kozek [3] found that for elliptical symmetry also the 
joint optimization (g and 7) results in Hermite functions 5 . 
Explicitely known is the joint optimum for C being Gaussian 
[5]. In (13) there is an invariance with respect to (affine) 
canonical transformations. See for example [6] for a review 
on the metaplectic representation. Essentially, this means that 
we can translate, rotate and shear C into some prototype shape 
being canonical equivalent but with further symmetries. 

Focusing now on localization operators Lc. 7 with C := 
Xu/\U\ and let be A max := A max (ic, 7 )- One can show the 
following result: 

Lemma 2 Let U canonical equivalent to square. Then: 

erf(^\U\/l) 4 , . „ . , -m |rrhll n - 
—-2 < sup A max <min(e « , \U\ ) (15) 

The upper bound is general and holds for any shape. 

The proof can be found in [7]. The upper bound can also 
be found in [5]. The lower bound can be calculated for other 
shapes (like a disc) as well. As a result for Problem 1 it follows 
from Lemma 1 that: 

\ 1/ max(t>,2) 

\U\ -erf(v^74) 2 J (16) 

C. Bounds for SINR 

On the other hand, from the previous results we are also 
able to obtain for |A| < 1 an upper bound on SINR as 
follows. Assume that {7 n }ngz 2 establishes a frame. Then we 
can compute 7 = O7 to obtain a tight frame {~f^} n ei, 2 
with £> 7 o = |A| _1 . For a tight frame holds equality in 
(12). It remains to find the optimal g which can be found 

5 Kozek considered 3 = 7. However one can show that for elliptical 
symmetry around the origin the optimum has also this property. 



for example as the maximizing (generalized) eigenvalue as 
shown in [4] (see also Section V-A). Let us call SINR* := 
max 9 SINR(g, 7 , A). The upper bound in (15) of Lemma 2 
gives then for \U\ < e that: 

SINR* < ((a 2 + lAr 1 )^ - l) 1 (17) 

Effectively, the upper bound (17) represents the assumption 
that we would be able to retain with 7^7° the localization 
upper bounded by the rhs of (15), thus ignoring the loss due 
to the BLT. This approach is valid to obtain an upper bound. 

As already discussed at the end of Section III-E, in general 
no quantitative estimates exists on the localization loss, i.e. 
it is unknown to what degree (17) is achievable. Assuming 
negligible localization loss ("no BLT") the lhs of (15) gives 
us also that: 

no BLT / \U\ \ 

SINR* > f (c + A ) 1 .. -1 (18) 
\ erf(0r|C/|/4) 4 J 

For uncritical A (or A° respectively) the mapping O can be 
well conditioned (which still depends on 7) such that (18) 
is an estimate of achievable values of SINR*. Our numerical 
results, shown in Section V-C, support these assumptions for 
|A| = 1/2. 

V. Localization Algorithms and Numerical 
Verification 

A well-known approach is to adapt the lattice A and the 
time -/frequency spreads at and 07 of the pulses to the scatter- 
ing function C of an ensemble of WSSUS channels where the 
device should operate. In practice various environments have 
to be classified by their spread in delay (Ct) and mobility 
(spread Cj in two-sided Doppler frequencies) such that a 
potential mobile device will support a group of different 
modes {g,j, A}. For A = diag(T, F) this means to fulfill 
essentially T/F « a t /af m Ct/Cf which is known as "pulse 
and grid scaling" [3], [25], [4]. However, up to scaling (and 
displacing) this approach does not further fix g and 7. From 
time-frequency uncertainty reasons a Gaussian or the "tighten" 
Gaussian (the IOTA approach) can be used for example. More 
detailed properties of the scattering function can be exploited 
in pulse design with the following methods directly based in 
the theory presented so far. 

A. "Mountain Climbing" 

The general SINR problem is known to be a quasi-convex, 
convex-constrained maximization problem [4], that is a global 
optimization problem. The same holds, separately, for the 
localization problem: 

sup X max (L c , g ) = sup A max (L£ ) (19) 
9 7 

as well, which can be formulated as a convex, convex- 
constrained maximization problem. But, a so called "mountain 
climbing" algorithm can be used to perform a local optimiza- 
tion (called "Gain optimization" in [4]). Essentially this means 
to start with an appropriate g^ ~> (for example a Gaussian). 



In nth iteration step 7^™-' is calculated as the maximizing 
eigenfunction of L c 9 (»-i) ■ From 7W the operator L^, (n) is 
constructed to compute g' n ' as its maximizing eigenfunction. 
Similarly, the SINR can also be optimized directly with the 
notion of generalized eigenvalues (more details in [4]). 

B. A Solvable Lower Bound 

In contrast to the iteration method presented in the last 
section, it is possible to solve a lower bound of localization 
problem in one step. A simple convexity argument shows that: 

(C,|A 37 | 2 ) > |(. g ,Q 7 )| 2 (20) 

where Q is an operator with spreading function C. The 
optimum of the rhs of (20) over g and 7 are the maximizing 
eigenfunctions of QQ* and Q*Q (on the level of matrices 
computable by the SVD). 

C. Numerical Verification 

In the following we will compare the performance on 
localization and SINR of the different methods presented 
above. We consider a pulse shaped OQAM/OFDM system as 
presented in III-C, i.e. A = diag(T,F) with |A| = 1/2. The 
scattering function is Xu/\U\ for U = [0, Td] x [— Bd,Bd] 
and \U\ fw 0.29 which is a highly interference dominated 
scenario (but still underspread). We have varied the ratio 
between time and frequency dispersion (different Td/(2Bjj)) 
and A is always adapted according to the grid scaling rule 
(as explained above). All computations are performed on a 
discrete time system of length L = 512 (for more details 
see [4]) and bandwidth W, such that the equivalent discrete 
variables are r d ^ := T d ■ W and B ( j? ] := B D /W ■ L. While 
keeping constant (r d (£)) + 1)(2 ■ Bp* + 1) = 150 the discrete 
ratio R := (t^ + \)/{2B^ + 1) has been varied according 
to the following table: 
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In summary we will show the following methods: properly 
scaled Gaussians (Gauss), {,9°,7°} from the method of Sec- 
tion III-E with properly scaled Gaussians as input (IOTA), 
maximizers of the rhs of (20) (SVD), algorithm of Section 
V-A (localg). {g°,7°} from the method of Section III-E 
with "localg" as input (localg-tight) and the iterated SINR 
algorithm explained due to limited space only in [4] (sinralg). 
The last algorithm is in principle equivalent to the eigenvalue 
optimization in "localg" for particular noise level (here a 2 = 
— 20dB), however on the level of generalized eigenvalues. 

We will compare them to the theoretical estimates: the rhs 
in (17) (upper) and the rhs of (18) under the assumption that 
there will be no localization loss due to the BLT (lower). 

Fig.l shows the corresponding localization (9) and Fig. 2 
the SINR. It can be seen that the bounds are suited to describe 
the averaged performance of a multicarrier system in doubly- 
dispersive channels based only the support parameter \U\. It is 
interesting that also lower bound gives a quite useful estimate 
on the performance which is due to the uncritical condition 
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Fig. 1. Localization (C, | A 97 1 2 ) for \U\ » 0.29, lattice |A| = 1/2 and 
OQAM signaling. R has been varied. 
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Fig. 2. SINR /or \U\ » 0.29 a? noue power cr 2 = -20rfB, lattice |A| = 
1/2 and OQAM signaling. R has been varied. 

number of O at lattice density | vV | 1 = 2. For OFDM systems 
operating at "more critical" lattices much more loss to the 
BLT has to be expected. The inner values of R correspond to 
"full" doubly-dispersive channels where no common eigen- 
structure can exists such that each decomposition is always 
an approximation. The boundary values in turn are single- 
dispersive channels having a joint exact eigenstructure, i.e. it 
is possible to completely suppress interference. 

VI. Conclusions 

We have introduced the theory of pulse shaping with focus 
on the WSS US-averaged SINR and considered the question 
of approximate eigenstructure of time-varying channels with 
compactly supported spreading. With increasing demand on 
bandwidth and efficiency the understanding of the fundamental 
limits in both directions will be important for future wireless 
communication systems. For Weyl-Heisenberg signaling, as a 
general description of OFDM and OFDM/OQAM communica- 
tion systems, we found that both are on the level of localization 
equivalent. We have shown that simple localization bounds 
can be used to obtain general estimates on the eigenpair 
approximation behavior E% and the SINR itself. With the 
latter, for example, we are able to show fundamental limits on 
achievable performance based on simple statistical properties 
of the time-varying environment. 
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